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SIMULATION OF TURBULENT WALL PRESSURE 

By 

Robert L. Ash^ 

SUMMARY 


A computer program has been developed to simulate the transient 
wall pressure field produced by a low speed fully turbulent boundary 
layer. The theoretical basis for the simulation has been discussed 
and preliminary results from a pressure simulation are presented. 

INTRODUCTION 

The purpose of this report is to present a method for simulating 
the unsteady pressure fluctuations produced by a low speed incompress- 
ible turbulent boundary layer. Attention has been restricted to flow 
along flat plates with no pressure gradient. The simulation has been 
designed for utilization in analysis of compliant wall candidates for 
the Langley Research Center drag reduction program (refs. 1 through 5), 
A one-dimensional simulation is presented here which is compatible 
with both finite difference and finite element structural programs. 
Two-dimensional simulations can be developed by randomly accessing 
the one-dimensional simulation, but that development is not pre- 
sented here. 

Properties of the turbulent simulation have been extracted 
primarily from three experimental papers. The very extensive early 
work of Bull (ref. 6) has been employed to model the frequency 
spectrum, the convection velocity, and the decay of pressure 
fluctuations. Burton's (ref. 7) measurements ha’^e been used to 
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justify the assumption that turbulent sublayer "burst" data (see 
ref. 8 for a survey of basic sublayer phenomena) can be used to 
infer spatial and temporal distributions of pressure fluctuations 
in the outer portion of the boundary layer. (This inference does 
not necessarily imply that the burst structures ultimately form 
large scale pressure fluctuations.) Burton's work has also been 
used to estimate the distribution of pressure fluctuation amplitudes 
occurring in a collection of events, Offen and Kline's (ref. 9) 
burst distribution measurements have been used to model the spatial 
and temporal distribution of the pressure fluctuations. 

In order to minimize the influence of numerical resolution on 
the simulation, a random stepping procedure has been employed in 
space and time. The magnitude of the pressure fluctuation at a 
particular point at any time is constructed from a collection of 
the randomly generated pressure fluctuation "events". Each of the 
events which contribute to the local pressure level has been con- 
structed to preserve the turbulence statistics measured in experi- 
ments. Probability distributions for time between events, distance 
between events, fluctuation frequency, and amplitude have been 
employed. Each probability distribution has been assumed statisti- 
cally independent from the otheic, and a conventional Monte-Carlo 
approach has been used. Some detail is provided in the following 
section. 


MONTE-CARLO METHOD FOR SIMUIrATING WALL PRESSURE 
FREQUENCY SPECTRUM 


Bull (ref. 6) has reported the wall pressure frequency spectrum 
shown in figure 1 for a low speed turbulent boundary layer. He has 
shown his measurements are well represented by the empirical relation: 
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where is the dynamic pressure, 6* is the displacement thick- 

ness, U is the free stream velocity, and ui is the radian 

oo 

frequency. 

By assuming equation (1) adequately describes the frequency 
distributions of the pressure fluctuations occurring in a large 
number of events, the probability of a particular frequency occurr- 
ing in a particular event can be calculated. That probability would 
simply be the normalized foi n of equation (1) or: 


P(u)) = (J)p(w)/ 


'/ 


(|)p(a)) 


do) 


( 2 ) 


The probability distribution cannot be used directly in a simulation. 
Rather, the cumulative probability distribution must be used. The 
cumulative probability is defined as the probability an event will 
occur with a frequency less than or equal to w, whereas probab- 
ility was the likelihood an event would occur between o) and 

A 

0 ) + dcj. Probability P(w) and cumulative probability P(u) are 
related by 


P(lo) 



P (w) du) 


(3) 


However, in this case, P(w) is zero for all values of u less 
than zero, leaving: 


P(oj) =/ P((o)d(jj 
•'o 


(4) 


Equation (1) can be employed in (2) and (4) to write: 


P(o)) = 1 - 0.5916e"^^ - 0.5443e"‘^^^ + 0.1359e“®^ 


(5) 


where X 


Dimensionless frequency is shown as a function of 


0l)6 * 

U_ 


cumulative probability in figure 2, 


Most digital computers have available efficient random number 
generation subroutines. Generally, they produce uniformly distributed 
random numbers over the interval between zero and unity. If values 
of frequency are assigned from figure 2 , using randomly generated 
numbers between zero and one, the frequencies will be distributed 
in agreement with Bull's (ref. 6) data. A brief explanation of the 
basis for this approach is given in reference 10. 

Graphical procedures are not acceptable for efficient computer 
utilization. Furthermore, because equation (5) is a transcendental 
equation for X a function of P, it cannot be used efficiently in 
calculating X for a given random number. Since equation (5) is 
itself an empirical curve fit, this investigation has employed an 

A 

approximate equation for X a function of P to streamline computer 
calculations. The equation; 
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0.2173 
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appears to satisfactorily represent equation (5), as shown in figure 
2. Detectable errors occur only on the upper end of the plot between 

A 

0.8 and unity, with a maximum error of about 10 percent at P = 0,9, 
That error is within the uncertainty of Bull's (ref, 6) data fit. 

The same procedure has been used to generate the other turbulent 
characteristics simulated here. The remaining discussion will only 
document the appropriate probability distributions and their repre- 
sentations in the computer simulation. 
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Amplitude of a Pressure Fluctuation Event 

Experimental measurements of the amplitude of a pressure fluc- 
tuation produced by a single fluctuation event is not currently 
possible because background contributions from other events are 
always present. Such an amplitude is required in the present 
simulation. The only data which is related (indirectly) to indi- 
vidual fluctuation amplitudes are the measurements of Burton (ref. 

7) for the threshold pressure fluctuation which appeared to be 
responsible for a burst event. His measurements suggest the dis- 
tribution of amplitudes is Gaussian. Since Bull's (ref. 6) measure- 
ments indicate the root mean square (rms) pressure fluctuation is 
given by P^jj^g ~ speeds of interest, a Gaussian distri- 

bution can be employed using a standard deviation of 3 t^ about a 
zero average pressure. The Gaussian distribution was simulated using 
conventional techniques (ref. 11). 

Time Interval and Spacing Between 
Pressure Fluctuation Events 

Offen and Kline (ref. 9) have studied the time interval between 
burst events in the wall region of a low speed turbulent boundary 
layer. They have found that the time interval can be scaled with 
outer flow variables and the probability of a new event varies 
with dimensionless time 9 given by 6 = U^t/6*, as shown in 
figure 3. That data can be represented by gamma distribution 
function: 

p(0) = 6 ^ e"®"'^/[r(a + 1) (7) 

where a = 2.2 and B = 16.4, as shown in figure 3. 

Equation (7) was integrated numerically to obtain the cumulative 
probability distribution P(0) shown in figure 4. That distribution 
has been approximated by: 

0 = 32.2 + 72 p2 + 0.63 tan ^ P. (8) 

Pq + 0.619 ® 2 0 

and is also shown in figure 4. 


The spatial distance between events can be gotten by assuming 
time and space are related through the friction velocity u^« 
Examination of Of fen and Kline's (ref. 9) data has suggested that 
the Ax and At for spacing between events are related by: 
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Consequently, random numbers can be used in equation (8) to generate 
X where 


X 


&* 


u 



( 10 ) 


is the spacing between fluctuation events. 

Convection Velocity and Decay Rate 

Data from Bull's experiments (ref. 6) have been used to model 
both convecc"'’ on velocity and decay rate. Although Bull reports the 
convection veloc.’ty varies with distance from the fluctuation source, 
a constant value of convection velocity, u^, given by 

u^ = 0.8 , (11) 

c ® 


has been used in the present simulation. The constant assumption was 
made because of difficulties in the numerical calculations, but may 
be justified from other experimental data [see Willmarth (ref. 12) 
for example] . 


Bull's experiments (ref. 6) indicate spatial decay scales with 

^ XU 

+ T 

X = — ^ . His data has been used to vary the amplitude of the 
pressure fluctuation in a particular event. If A is the original 
amplitude, A(x ) has been modeled by; 


A(x^) 
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Both convection velocity and decay rate have been "modeled" 
rather than "simulated" in the sense that all events are assumed 
to have the same convection and decay properties. At this point 
equations have been developed to simulate frequency, amplitude, 
and spacing of pressure fluctuation events. The remaining dis- 
cussion is intended to explain how these models have been put 
together to simulate a turbulent wall pressure. 

SIMULATION OF THE TURBULENT WALL PRESSURE 

The ultimate output of this simulation is to be a simulated’ 
pressure field over a specific distance (test model) with a specific 
time interval. In order to generate that simulation, both a develop- 
ment length and a start-up time are required. Theoretically, dis- 
turbances extending upstream to infinity can contribute to local 
pressure fluctuations which implies that the start-up time would 
also be infinite. Practically, a finite start-up length can be 
used which accounts for nearly all of the local pressure fluctu- 
ation. However, because of the random time stepping procedure 
employed here, the start-up time had to be determined by numerical 
testing. 

The decay rate defined in equation (12) can be used to specify 
the development length. If fluctuations less than one percent of 
their original value are neglected, then the start-up length is 
given by: 

xu 

X = 425,000 = — - (13) 

V 


If u =1.5 m/sec and v = 1.5 x 10“^ m^/sec, the start-up length 
T 

is 4.25 m. 

The procedure employed here has been to define the origin (x = 0) 
as the front of the start-up length. Assuming the start-up length is 
x_, the front of the model would be at x = x^. Using random num- 
bers, frequency f, amplitude P^, origin Ax, and time At are 
generated for a particular disturbance using the previously described 
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equationf,. The amplitude of the pressure fluctuation can be positive 
or negative; all other quantities are positive. 

A single-cycle, sinusoidal pressure fluctuation represents a 
single event. If is the location of the previous pressure 

fluctuation, the new fluctuation is located at 


X = X + Ax 

p 


(14) 


The time over which the event is sensed at a particular spatial loca- 
tion T is given by: 


rp _ 1 _ 27T 

f ~ U) 


(15) 


Using the ronvection velocity u^ the spatial length of the event 
is given by: 


= u T 
D c 


(16) 


and the origin of the disturbance is assumed to occur at 

X =x-W_=x +Ax-uT , (17) 

o D p c 

The time of birth of the event, t^^, is assumed given by 

tfa = + At (18) 

where t is a reference time which will be discussed later, 
o 

If the disturbance is assumed "born" instantly over the interval 

between x and x, at time t, , the distance traveled by the dis- 
o b 

turbance at some later time t is given by: 

D = u^ • (t - t. ) (19) 

c b 
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and the decayed amplitude P (t) is: 

d 


(t) 

a 
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( 20 ) 


If Xg is a particular location on the test surface where the 
pressure history is desired, then the arrival time for the disturbance 
is : 


t = 
a 


+ 


X - X 

s 


u 

c 


and the departure time is: 


t 


d 



( 21 ) 


( 22 ) 


Furthermore, if t is any time in the interval, t < t < t. , the 

a ■■ “ b 

pressure produced at x^ by that particular event is: 

t - t 

P'(x^, t) = P^(t) sin 2tt (23) 

S a 


Equation (23) then represents the pressure contribution at desired 

location x at time t and is the basis for all further calcula- 
s 

tions. 

Assuming reference x... t^ is )cnown, pressure fluctuations 

of the type just described can be generated at a random set of 

points, starting at x = 0 and stepping to x = x_ + x , where 

D m 

x^ is the model length, using randomly generated Ax's. In 
addition, times of birth can be distributed randomly above t^ 
by using random At's in equation (18). Then, using equation (23) 
the pressure contributions at particular points can be calculated. 

If At is defined as the average time step generated in 

d 

stepping from x=0 to x=x_+x, when the end of the model 

D m 

is reached, calculations can return to the origin and step down 
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the model again. The reference time would then be updated a&: 


t 

o 


t 

o 


+ At 

a 


( 24 ) 


That process can be repeated over and over until a particular 
time is reached (t = t ) . Assuming that part of the simulation 
is completed, two problems still remain — how to specify to and how 
to store the simulation. 


It is reasonable to assume the pressure simulation should be 


over the time interval 


0 < t < t 
— - max 


However, t cannot be set 


equal to zero to start with, because start-up transients would be 
present in the simulation. Since the time required for a disturb- 


ance to travel from x = 0 
start-up time would be: 


to 


X = X_ + X 

D m 


IS 


X_ + X 

D m 
u 


a logical 


t . 
min 


x_ + X 
D m 

u 


(25) 


However, t . should be larger because some of the disturbances may 
min ^ ^ 

occur upstream from x = 0, i.e., x^ can be negative. A conser- 

vative value for start-up time given by 


t 

o 


-1.44 


t . 
min 


(26) 


has been used in the present simulation. Then at t = 0, start-up 
transients should no longer exist. 


Locations and time steps used in the storage of the simulated 

pressure have no effect on the calculations and are specified based 

on the spatial and time resolutions for the structural calculations. 

An Xg array can be specified, as well as the desired structurally 

compatible time step ^^ 5 * Then, if a single disturbance passes over 

X between t and t,, integer start number and stop number 

s a d ^ go 

N ^ given by: 

stop ^ ^ 
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and 


(27) 




0.99 


N 


Stop 



can be used to assign storage time locations for P' t) . That 

is, the pressure P(x , t ) has been increased by an amount 

s n 

P' (x , t ) for t = nAt , N < n < N . . It is important to 

s' n n s go " - stop ^ 

note that P'(x , t ) is not the pressure because more than one 
s n 

event may be over at a given time. 

A program for performing the pressure simulation just described 
is included in Appendix A. Output from a sample simulation run is 
shown in figure 5. It shows the variation of pressure with time at 
a single point. All aspects of the simulation have not been verified 
at this time. That is, output power spectra and correlations have 
not been examined and compared with experimental data. However, the 
root mean square pressure fluctuation and average pressure agree 
quite well w.*-n their desired values (3 t and 0, respectively). 
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FORTRAN PROGRAM FOR SIMULATING 
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